A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

04(11.09.1 Dwl, 11.09.1 Dw2, 11.09.03, 11.11.1, 11.12.1, 11.19.2) 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



> 
o 



Neutral hydrogen in the nearby galaxies Dwingeloo 1 and 
Dwingeloo 2 

W.B. Burton^, M.A.W. Verheijen^, R.C. Kraan Korteweg^ and P.A. Henning'' 

^ Sterrewaclit Leiden, Postbus 9513, 2300 RA, Leiden, The Netherlands 

Kapteyn Astronomical Institute, Postbus 800, 9700 AV, Groningen, The Netherlands 
Observatoire de Paris-Meudon, DAEC, 5 Place Jules Janssen, F-92195 Meudon-Cedex, France 
Department of Physics and Astronomy, University of New Mexico, Albuquerque, NM 87801, U.S.A. 



Received . . . , 1995; accepted 



1995 



> 

o 

O 



in 



Oh; 

6 ■ 



X 



Abstract. We present observations made with the West- 
erbork Synthesis Radio Telescope of H i emission from 
Dwingeloo 1, a nearby barred spiral discovered during the 
Dwingeloo Obscured Galaxies Survey for galaxies hidden 
in the Zone of Avoidance, and of Dwingeloo 2, a small 
galaxy discovered in the beam of these WSRT observa- 
tions. The WSRT data reveal the position of the dynam- 
ical center of Dwl, its systemic LSR velocity, its total 
(projected) width in velocity, its inclination on the sky, 
and its integrated Hi flux, as well as details of the veloc- 
ity field and gas distribution. Dwl is the nearest grand- 
design barred spiral system, and is probably amongst the 
ten largest galaxies closer than about 5 Mpc. We report 
here also the discovery of Dwingeloo 2, a small galaxy 
located within the WSRT primary beam as pointed to 
Dwl. In view of its angular and kinematic proximity to 
Dwl, Dw2 may well be a companion to the larger system. 
The two galaxies are probably both members of the group 
containing Maffei 1 & 2 and IC342 and may influence the 
peculiar motions within that group and the morphology 
of its individual members. 

Key words: galaxies: individual: Dwl - galaxies: individ- 
ual: Dw2 - galaxies: interstellar medium - galaxies: spi- 
ral - galaxies: kinematics and dynamics - galaxies: Local 
Group 



1. Introduction; the Dwingeloo Obscured Galaixies 
Survey 

Some 20% of the optical sky is significantly obscured by 
absorption and scattering of light by intervening interstel- 
lar material in the Milky Way. As a consequence of this 
obscuration, the inventory of nearby galaxies is certainly 
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incomplete closer than about 5° to the galactic equator, 
where the optical and infrared sky is grossly confused 
by local stars and other objects and where deep optical 
searches for obscured galaxies fail (see Kraan-Korteweg 
& Woudt, 1994). The boundaries of the Zone of Avoid- 
ance are difficult to quantify; there are patches of sky 
beyond |6| ~ 5° where the optical depth is also substan- 
tial. The intervening interstellar medium which is largely 
opaque at optical wavelengths is transparent at the 21-cm 
wavelength of neutral hydrogen. Reviews of the use of the 
H I line to search for galaxies in the Zone of Avoidance 
have been given by Kerr and Henning (1987), by Kerr 
(1994), and by Henning (1992, 1994). 

The Dwingeloo Obscured Galaxies Survey has been 
searching for H i emission from galaxies hidden behind the 
obscuring material causing the Zone of Avoidance, using 
the 25-m radio telescope of the Netherlands Foundation 
for Research in Astronomy, located in Dwingeloo. The 
project is motivated by the hope to find the H i signatures 
of fairly nearby galaxies in the velocity range — 4000 km 
s~^. Such detections would contribute to investigations of 
the dynamics of the Local Group, of the morphological 
distribution of nearby galaxies, of the age of the Universe 
as determined by the timing argument, and of the den- 
sity parameter fl as revealed by the velocity of the Local 
Group relative to the Microwave Background. 

It is expected that the Dwingeloo Obscured Galax- 
ies Survey will continue, utilizing the telescope essentially 
full-time, for a duration of several years. Before this search 
began, the Dwingeloo 25-m telescope had been used over 
a five-year period to produce the Leiden/Dwingeloo at- 
las of H I in our Galaxy over the entire sky at ^ > —30° 
and within the velocity range —450 < v < +400 km s^^ 
(Hartmann, 1994; Hartmann & Burton, 1995). That sur- 
vey represented the first use of the 25-m telescope with the 
new DAS 1000-channel digital autocorrelator spectrome- 
ter combined with the receiver, characterized by a system 
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Fig. 1. Mosaic of channel maps showing the Hi emission from Dwingeloo 1, after subtraction of the 21-cm continuum radiation 
and after application of the CLEAN algorithm, but uncorrected for primary-beam attenuation. For each panel the appropriate 
central LSR velocity is indicated; the cross indicates the position of the dynamical center. The ellipse in the upper-lefthand 
panel indicates the approximate size and orientation of the optical image obtained with the Isaac Newton Telescope in the 
/-band image (see Loan et al. 1995). Contour lines in the channel maps correspond to —0.42, —0.21, 0.21 (the 2-cr level), 0.42, 
0.85, 1.3, 1.7, 2.5, 3.4, and 4.2 xlO'^" Hi atoms cm~^. The first five channel maps show contaminating Hi emission from our 
own Galaxy. The lower-righthand panel shows the {\m-CLEANeA) continuum radiation from the field, for which the individual 
channel maps have been corrected. Contour lines in this continuum map correspond to —1.3, —0.62, 0.62 (the 2-cr level), 1.3, 
2.5, 5.0, 7.5, 10, 13, and 15 K; dashed-line contours represent the negative values. 



temperature of about 35 K, which is a prototype of the 
sort mounted on the telescopes of the Westerbork array. 
The experience gained with the Leiden/Dwingeloo Milky 
Way survey confirmed that the interference environment 
was benign enough, and the telescope electronics stable 
and sensitive enough, to justify the search for obscured, 
relatively nearby, galaxies. 

The performance and operating configuration of the 
Dwingeloo 25-m telescope currently pertaining is that de- 
scribed by Hartmann (1994), Burton & Hartmann (1994), 
and Hartmann & Burton (1995). For the obscured-galaxies 
project, the radio receiver is tuned to a central frequency 
corresponding to a velocity of -1-2000 km s~^ above the 
rest frequency of the H I hyperfine transition; the band- 
width of the 1000-channcl spectrometer is set to cover the 
range Q < v < -1-4000 km s^^ at 4 km s~^ resolution. (Ra- 
dial velocities are defined here following the radio conven- 
tion and arc referred to the Local Standard of Rest de- 
fined by the Standard Solar Motion of 20 km s~^ towards 
a,5 = 18'', -1-30°, epoch 1900.) 

The Dwingeloo Obscured Galaxies Survey involves a 
strategy with several different facets. Two blind searches 
will fully sample the Zone of Avoidance; first, a program 
of short (5-minute) integrations has been searching for 
H I emission from large, or particularly nearby, hidden 
galaxies; a program of longer (12 x 5-minutes) observa- 
tions has now been initiated looking for smaller, or more 
distant systems. Both blind searches involve observing on 
a honeycomb grid with lattice points separated by 0?35; 
because the FWHM beam of the 25-m telescope operating 
at a wavelength of 21 cm is 36 arcminutes, the honeycomb 
grid is fully sampling the beam. The blind searches are 
covering the regions of full optical obscuration at \b\ < 5°, 
over most of the accessible range in longitude. The initial 
period of the survey has been directed toward the longi- 
tude range in the second quadrant where the Supcrgalactic 
equator crosses the Zone of Avoidance and where a large 
population of galaxies is expected a priori. (Most of the 
galaxies all of them previously known ones identified as 
contaminating nuisances in the Leiden/Dwingeloo survey 
of H I in the Milky Way occur in a vertical swath crossing 
the equator of the Milky Way near the center of the sec- 
ond longitude quadrant: see e.g. figure 4.4 in Hartmann, 
1994, or figure 5 in Burton & Hartmann, 1994.) Although 
the blind searches are the principal facets of the strategy, 



a limited directed search is being carried out involving 
longer integrations toward candidate objects selected from 
optical and infrared datasets. Integration times of an hour 
per spectrum yield sensitivities of about 0.01 K over the 
range < < -1-4000 km s~^. 

Many problems can plague single-dish 21-cm spectra. 
In particular, the expected emission signature of an ex- 
ternal galaxy can be mimicked by interference signals of 
unusual type, i.e. not appearing as a single-channel spike 
or as a sin(a:)/a; ringing (see figure 2.13 of Hartmann, 
1994); the signal from a low- velocity galaxy like either 
Dwl or Dw2, or like Maffei 2, might be at least partly 
blended with that from conventionally-behaving gas as- 
sociated with the Milky Way, or it might resemble that 
of a (very) high-velocity cloud associated with the Milky 
Way; the signal might be confused with edge-of-bandpass 
effects; or the signal might be contaminated by enhanced 
system temperatures due to the appearance of a Milky 
Way H II region, a plethora of which occur at the latitudes 
being searched (see e.g. Lockman's 1989 catalog). Each of 
these problems could be dealt with at the 25-meter tele- 
scope itself, but an important aspect of the strategy of the 
Dwingeloo Obscured Galaxies Survey is the intention to 
seek rapid independent confirmation of suspected detec- 
tions by using the Westerbork Synthesis Radio Telescope 
in the "snapshot" mode. Thus the detection spectrum of 
Dwingeloo 1 (sec Kraan Korteweg et al., 1994) was dis- 
cussed amongst members of the project team in a meeting 
in Dwingeloo on Thursday, August 4, 1994, and by the end 
of the day on Monday, August 8, a WSRT snapshot was 
in hand confirming the detection and giving preliminary 
information on the H I position and flux. By the end of 
that week, optical and near- as well as far-infrared iden- 
tifications of the newly detected galaxy had also been ob- 
tained (see Kraan-Korteweg et al. 1994; Loan et al. 1995). 
We discuss here full 4 x 12 hour synthesis observations of 
Dwingeloo 1, and the discovery in these observations of 
another galaxy, which we call Dwingeloo 2 and which is 
probably a companion to Dwl. 

2. WSRT observing program 

The Westerbork Synthesis Radio Telescope is an interfer- 
ometer comprising an array of 14 telescopes, each with a 
diameter of 25 meters. The east-west alignment implies 
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that the technique of earth-rotation aperture synthesis re- 
quires an observation of 12 hours duration to obtain a syn- 
thesized field. The observations reported here were carried 
out in five sessions during August and September, 1994, 
and resulted in a full synthesis of 4.4 x 12 hours duration, 
including the snapshot material. The four eastern, mov- 
able telescopes were set at different positions during the 
different 12-hour periods, and the signals from the fixed- 
movable telescope pairs were correlated; this provided a 
total of 152 different baselines ranging in length from 36 
to 2754 meters with a standard increment of 18 meters. 
As a result, the instrumental grating rings caused by the 
constant increment in baseline length lie well beyond the 
FWHM of the primary beam and do not interfere with 
extended emission. 

The synthesized beam has an elliptical shape with di- 
mensions a X S = 11.8 X 13.8 arcsec. The FWHM size 
of the primary beam of the WSRT when operated at a 
wavelength of 21 cm is 37^6, quite similar to that of the 
single-dish Dwingeloo telescope. A series of spectra made 
with the 25-m telescope at neighboring pointings around 
the source later to be identified as Dwingeloo 1 had given 
a rough indication of the direction towards which the 
WSRT should be pointed for the initial "snapshot" obser- 
vation, and the suspicion (correct, as it turned out) that 
the smudge on the red Palomar Sky Survey plate (^^936) 
catalogued by Hau et al. (1995) might be the counter- 
part to the H I feature, led to centering the WSRT on the 
direction a ^ 02''53™02^ S ^ +58°42'36" (epoch 1950). 
The central frequencies of the WSRT receivers were tuned 
to correspond to an LSR velocity of 106 km s~^. The 
single-dish detection had shown the double-horned signa- 
ture characteristic of H i emission from a spiral galaxy, 
but the extent to which the lower-velocity horn was con- 
taminated by emission from the Milky Way was not clear 
until the WSRT data confirmed that the central-frequency 
tuning of the receivers was reasonably accurate. A total 
bandwidth of 2.5 MHz was covered by 63 channels; during 
the observations an on-line uniform frequency taper was 
applied, yielding a velocity resolution of 9.89 km s~^. 

The interferometric [/T^-plane data were calibrated, 
fiagged if necessary, and Fourier transformed with the 
NEWSTAR software package. Before Fourier transform- 
ing, the J/V-plane data were convolved with an expsinc 
function to reduce aliasing effects, and a Gaussian base- 
line taper was applied to suppress the sidelobes of the syn- 
thesized beam. Furthermore, the UV-data were weighted 
according to the local density of data points in the UV- 
plane. Due to decreasing sensitivity near the edges of the 
bandpass, the first 3 as well as the last 7 channels were dis- 
carded. This procedure resulted in a data cube containing 
1024 X 1024 X 53 pixels, covering a field of angular di- 
mensions 1?42 X 1?66 with a synthesized beam of FWHM 
of 11'.'8 in a and 13'.'8 in 6. The total velocity coverage 
was 428.7 km s~^. Application of the various tapers and 
convolution functions yielded an rms noise in the channel 



Table 1. Westerbork observing parameters for the Dwingeloo 
1 field 



field center (1950.0) a 




02''53™02'' 


5 




58°42'36" 


dates of observations 




8 & 12 Aug. 1994 






6, 26 & 29 Sept. 1994 


length of observations 




4.4x12 hours 


number of different baselines 




152 


baselines (min-max-incr) 




36-2754-18 meters 


FWHM of synthesized beam (a 


X 5) 


11.8 X 13.8 arcsec 


radius of first grating ring (a x 


5) 


40.1 X 47.0 arcmin 


FWHM of primary beam 




37.6 arcmin 


rms noise in one channel map 




4.8 K 


central frequency 




1420.01 MHz 


central velocity {LSR) 




106 km s"^ 


bandwidth 




2.5 MHz 


number of channels 




63 


channel separation 




8.24 km s"^ 


velocity resolution 




9.89 km s"^ 


K=>mJy conversion, 






equivalent of 1 mjy/beam 




4.24 K 



maps of 1.2 mJy per beam. A set of 5 maps with the an- 
tenna patterns at different frequencies was constructed for 
use when applying the CLEAN algorithm to the channel 
maps. 

Further reduction and analysis of the data utilized the 
Groningen Image Processing SYstem (GLPSY), following 
practices described by Verheijen (1996). First, the data 
cube was smoothed to spatial resolutions of 30 and 60 arc- 
seconds in order to obtain higher sensitivity to extended 
emission at low column densities. These smoothed data 
cubes showed that H I emission from Dwl covers the ve- 
locity interval extending from 7 to 221 km s~^; the emis- 
sion from Dw2 extends from 32 to 156 km s~^. Because 
the channels at velocities below about 25 km s^-'^ are con- 
taminated by H I emission from the Milky Way, a map of 
the 21-cm continuum emission was made by averaging 10 
channels at velocities from 229.7 to 303.9 km s"-'^, which 
were judged to be free of any Hi emission either from the 
Milky Way or from the galaxies in question. The rms noise 
in the resulting high-resolution contimmm map is 0.37 
mJy/beam. The continuum map was subtracted from each 
channel map in the data cube. The continuum-subtracted 
channel maps which do contain H i emission from Dwl and 
Dw2, as well as the continuum map itself, were CLEANed 
by applying that algorithm in the standard manner (see 
Hogbom 1974). The map areas to be cleaned were deter- 
mined by the 2-cr contour in the un-CLi?^iVed map at 
60" X 60" resolution and extended by 1 arcmin in order to 
accoimt for possible emission suppressed by the sidelobes 
below the formal 2-a level. Note that these areas vary 
from channel to channel due to the rotation of both galax- 
ies. The selected areas were CLEANed down to the 0.3-cr 
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Fig. 2. Global kinematic and spatial properties of the Hi distribution in Dwingeloo 1. Panel a) shows, for comparison with 
the global Hi situation, a grcyscalc representation of the INT 7-band image (see Loan et al. 1995). The size and orientation of 
the optical image is represented by the ellipse in panel b), which shows the CLEANed continuum map, with contours drawn at 
levels of —1.2, —0.58, 0.58 (the 2-a level), 1.2, 2.3, 4.6, 9.3, and 19 K. Panel c) shows the global Hi line profile in mjy, corrected 
for the sensitivity of the primary beam; the vertical arrow shows the systemic velocity of Dwl, 110.3 km , as derived from 
the velocity-field data. Panel d) shows the deprojected Hi surface density in units of Mg pc^'^, as determined in tilted, rotating 
rings separately for the approaching and receding sides of Dwl; the full-drawn line corresponds to the average of the surface 
densities measured on the two sides; the vertical arrow shows the angular extent of the semi-major axis of the optical image. 
The panels labelled e) show position- velocity maps of the H i emission observed, at the indicated position angles, along the 
kinematic major axis (left) and along the kinematic minor axis (right) of Dwl. Contours in the major-axis position- velocity 
map correspond to -0.36, -0.18 (the 2-a level), 0.18, 0.36, 0.73, 1.1, 1.5, 2.2, 2.9, and 3.6 xlO^° Hi atoms cm"^; in the 
minor-axis map, the contours correspond to -0.42, -0.21, 0.21 {2a), 0.42, 0.85, 1.3, 1.7, 2.5, and 3.4 xlO^° Hi atoms cm"^. 
The dots follow the rotation curve as projected onto this position-velocity maps. In both panels, the horizontal dashed line gives 
the systemic velocity while the vertical one corresponds to the location of the center of rotation; the Milky Way contributes the 
emission contaminating the lower-velocity portion of each of these panels. The cross in the lower-lefthand corner of each of the 
position-velocity slices gives the spatial and kinematic resolution of the WSRT data. Panel g) shows the observed velocity field 
determined for Dwl by Gaussian fitting to individual H i spectra. The lighter greyscale and the black contours correspond to 
the approaching side; the darker greyscale and the white contours, to the receding. The heavily-drawn black contour follows the 
systemic velocity. The interval between contour lines is 20 km s~^, measured from the systemic velocity of 110.3 km s~^. Panel 
i) shows the modelled velocity field as a "spider diagram" derived by fitting the observed field with tilted rings in which the 
Hi describes circular rotation; the isovelocity contours in this panel are drawn as in g). Panel j) shows the difference between 
the observed and modelled velocity fields; the contours give the velocity residuals at levels of —10, —5, 5, and 10 km s~^. Panel 
h) shows the distribution of the total Hi column density across Dwl. As indicated in the text, the noise level varies across 
this map. Greyscaled pixels in the total-Hl map all have a signal-to-noise ratio > 1.5; the contours correspond to 1.8 (mean 
S/N = 3), 3.0, 6.1, 9.1, 12, and 18 xlO'^" Hi atoms cm^'^. The three panels in /) show the radial variations of, respectively 
beginning with the upper of the three, the inclination of tilted rings fit to Dwl, the position angle of these rings, and the 
rotation velocity. The rotation velocity was determined separately for the approaching and receding sides, as indicated. The 
full-drawn line in the rotation-curve panel does not correspond to the average of these two solutions, but to a separate solution 
incorporating complete rings. The filled circle in the two upper panels and in the one on the lower left indicates the FWHM of 
the synthesized WSRT beam. 



level. The CL£'^ A/' components were restored by Gaussian 
beams of FWHM dimensions 11" x 13", 30" x 30", and 
60" X 60", respectively. The rms noise in the continuum- 
subtracted and CLEANed channel maps is 1.2, 1.6, and 
2.3 mjy/beam for the maps at high, medium, and low 
resolutions, respectively. After CLEAMng, the map areas 
containing H I emission were redefined on the basis of the 
2-a level in the CLEANed maps at 60" x 60" resolution. 
The maps at 30" x 30" resolution give the best compro- 
mise between resolution and sensitivity and are thus the 
ones considered further here. 

The global Hi profiles of Dwl and Dw2 were deter- 
mined by summation of the total signal, primary-beam- 
corrcctcd, in those areas. Taking the primary beam atten- 
uation into account, the sum of the two global H i profiles, 
one for Dwl and one for Dw2, corresponds very well to 
the single-dish detection spectrum measured toward Dwl 
on the Dwingeloo 25-m telescope by Kraan-Korteweg et 
al. (1994), as well as with the spectrum observed subse- 
quently by Huchtmeier et al. (1995) using the Effelsberg 
100-m telescope . 

Relevant parameters of the WSRT observing program 
are summarized in Table 1. 



3. H I in Dwingeloo 1 

Fig. 1 is a mosaic of channel maps showing the H I emission 

from Dwingeloo 1; the panels represent alternate channels, 
each corresponding to the emission over an interval 8.2 km 
s~^ wide centered at the indicated velocity. The first five 
panels in the scries show contamination from H i emission 
from the Milky Way. The last panel shows the subtracted 
continuum map. H I emission from Dwl is evident in each 
of the displayed panels representing channels centered be- 
tween 15 and 213 km s~^. The emission shows the "spider- 
diagram" isovelocity characteristics which are the signa- 
ture of H I in an inclined spiral galaxy. The ellipse drawn 
in the upper-lefthand panel of Fig. 1 represents schemati- 
cally the orientation and size of the optical image of Dwl 
observed in the /-band with the Isaac Newton Telescope 
on La Palma (see Kraan-Korteweg et al. 1995; Loan et al. 
1995); the image itself is reproduced in Fig. 2a. The spatial 
extent of the Dwl H i channel-map slices is clearly larger 
than the extent of the optical-image ellipse. Although such 
a situation is common for spiral galaxies, we note that the 
optical obscuration at 6 = —0? 11 is so severe that the outer 
boundary of the /-band image is probably not a quantita- 
tive measure of the actual extended stellar dimensions of 
the system. 
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Fig. 3. Mosaic of channel maps showing the Hi emission from Dwingeloo 2, after subtraction of the background continuum 
radiation and after application of the CLEAN algorithm, but uncorrected for primary-beam attenuation. For each channel 
map the appropriate central LSR velocity is separately indicated; a cross indicates the position of the dynamical center. The 
lower- righthand panel shows the {un- CLEANed) continuum radiation from the field, which was subtracted from the individual 
channel maps. Contour lines in the channel maps are drawn as in Figure 1. Contour lines in the continuum map correspond to 
— 1.3, —0.62, 0.62 (the 2-cr level), 1.3, 2.5, 5.0, 7.5, 10, 13, and 15 K; the dashed-line contours represent the negative values. 



Fig. 2 shows the principal large-scale kinematic and 
spatial properties of the neutral hydrogen distribution 
in Dwingeloo 1. The global Hi profile shown in Fig. 2c 
was obtained by integrating the individual continuum- 
subtracted channel maps, giving the Hi flux density as 
a function of radial velocity. The global profile shows the 
two-horned signature characteristic of an inclined spiral 
galaxy, with the flux enhanced at the edgew of the pro- 
file by velocity crowding. The arrow in the global-profile 
panel indicates the systemic velocity of 110.3 ± 0.4 km 
s~^ as derived from the velocity field of Dwl. The radial 
velocity of Dwl with respect to the center of our Galaxy 
is Vg = 110 -|- 220cos(6) sin(/) = 257 km s~^, assuming a 
circular velocity of the LSR of 220 km s~^. 

A map of the spatial distribution of the total integrated 
Hi column density in Dwl is shown in Fig. 2h. This total- 
Hi map was constructed by adding all the channel maps, 
after having set the pixels outside the defined areas to 
zero. Because the number of added channels thus varies 
from pixel to pixel, the noise in the total-H I map varies 
as well from pixel to pixel, independent of the strength of 
the signal. 

The panels of Fig. 2e show position- velocity diagrams 

corresponding, respectively, to the major axis of the total- 
Hi map, determined at a position angle of 112°, and to 
the minor axis of the total-H I map, at a position angle of 
202°. 

The observed (projected) velocity field of Dwl shown 
in Fig. 2g was determined by fitting single Gaussians to 
each velocity profile through the data cube. Only those 
Gaussian solutions were accepted which met the three- 
fold criterion of returning a peak inside the area defining 
the H I emission, of returning a peak amplitude greater 
than 3 times the noise in the profile in question, and of 
yielding a peak velocity with an error less than a third of 
the velocity resolution of the WSRT data. 

The kinematic and spatial characteristics of Dwl men- 
tioned above are indicative of a well-behaved system. Such 
behavior justifies further analysis based on a model of 
tilted, rotating rings. This analysis technique has been 
described by Begeman (1989). Fig. 2i shows the results 
of modelling the observed velocity field with such tilted 
rings. The width of each ring was held constant at value 
equal to 2/3 of the size of the major axis of the synthesized 
WSRT beam. 

The ring fitting involved three steps. First, the sys- 
temic velocity, the center of rotation, and the rotation 
velocity were modelled for each ring keeping the inclina- 



tion and position angle fixed at values derived from fitting 
ellipses to the total H i map. No significant trend in either 

systemic velocity or in rotation center as a function of ra- 
dius was detected; in the subsequent steps, the rings were 
therefore treated as spatially and kinematically concen- 
tric. Second, the rotation velocity, the position angle, and 
the inclination were fit for each ring keeping the systemic 
velocity and center fixed at the previously-determined av- 
erage values. No significant trend of inclination from the 
mean value of 51° (measured such that an edge-on galaxy 
is inclined 90°) was detected as a function of radius which 
might be taken as indication of a morphological property 
of the galaxy such as warping; some localized wiggles in 
the parameters might be due to the influence of spiral 
structure. A slight change of position angle detected near 
the center of Dwl may well be caused by non-circular mo- 
tion due to the presence of the bar evident on infrared and 
optical images of the galaxy (see Kraan -Korteweg et al. 
1994; Loan et al. 1995). Finally, the rotation velocity and 
the position angle were fit for each ring, keeping all other 
parameters fixed. This last step was done in half rings, 
separately for the receding and for the approaching sides 
of Dwl, in order to investigate the degree of symmetry of 
the rotating disk. The substantial symmetry of the sys- 
tem is indicated by the two panels of Fig. 2e, where the 
mean rotation curve is projected on the position-velocity 
diagrams. 

The orientations of the rings as derived by the tilted- 
ring model were used to derive a radial H i profile by aver- 
aging the H I column density from the total H I map in an- 
nuli corresponding to the projected rings. The width of the 
annuli was set to 2/3 of the major axis of the synthesized 
WSRT beam. Fig. 2d shows the variation of the depro- 
jected H i surface density in solar-mass units with angular 
distance from the galactic center, plotted separately for 
the receding and the approaching sides of Dwl. The ver- 
tical arrow indicates the angular extent of the semi-major 
axis of the /-band image. 

Finally, Fig. 2b shows the 21-cm continuum emission 
from Dwl. This emission, which is contributed from the 
entire extended disk, has an integrated flux of 48.2 ± 7.5 
mJy and suggests on-going star formation in the galaxy. 

4. H I in the newly discovered galaxy Dwingeloo 2 

The Westerbork observations resulted in the discovery of 
a small galaxy, located some 21 arcminutes to the west- 
northwest of Dwl and quite likely its companion. The 
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Hi signal from this galaxy, which we call Dwingeloo 2, 
would have been heavily blended in the single-dish detec- 
tion spectrum of Dwl. The signals from the two systems 
are centered less than 20 km s~^ apart in velocity, and 
their angular separation is less than the FWHM of the 
beam of the Dwingeloo 25-m telescope. 

We were unable to find Dwingeloo 2 on the POSS I 
plates, but inspection of the red POSS II plate revealed a 
feature of low surface brightness, elongated like the WSRT 
Hi image. C. Aspin, and R. Tilanus have detected the 
stellar component of Dw2 in observations made using the 
IRCAM3 (see Puxley et al. 1994) camera operating at a 
wavelength of 2.2/im on the UKIRT facility on Mauna 
Kea. They have kindly provided the i^T-band image repro- 
duced in Fig. 4a. This image reaches a depth of rriK = 20; 
although the signal from the stellar component of Dw2 is 
weak, it has evidently been detected as a smudge inclined 
on the sky in the same manner as the gaseous component 
and as the POSS H feature. (We note that the possibili- 
ties now becoming available of deep imaging at 2.2/im, as 
well as of surveying at this wavelength, will likely lead to 
upward revisions of the sizes and stellar masses of heavily 
obscured galaxies, as well as to the discovery of additional 
ones.) 

Fig. 3 is a mosaic showing the H l emission from Dw2 

in individual channel maps, each centered at the indicated 
velocity. Despite the decreased sensitivity (40% of maxi- 
mum) of the primary beam at the off-axis position of Dw2, 
the small galaxy contributes easily-detected signals to the 
panels between 32 and 155 km s~^. The "spider-diagram" 
signature so pronounced in Dwl does not characterize the 
Dw2 channel maps: the H i disk of Dw2 apparently rotates 
as a solid body. We note, however, that the H i distribution 
in each channel map is not quite perpendicular to the 
overall kinematic major axis as would be expected for an 
unperturbed disk rotating as a solid body. 

Fig. 4 shows the principal kinematic and spatial global 
properties of the Hi distribution in Dw2; these properties 
were extracted from the WSRT data in the same manner 
as described for Dwl. The global Hi profile shown in Fig. 
4c docs not show the two-horned signature, because of the 
absence of the velocity-crowding which would follow from 
an extensive region spanned by an essentially flat rotation 
curve. The systemic velocity of Dw2 as derived from the 
velocity field is 94.0 ± 1.5 km s~^ and is indicated by the 
arrow in the global-profile panel. The radial velocity of 
Dw2 with respect to the center of our Galaxy is vg = 
94 + 220 cos(6) sin(/) = 241 km s'^ 

A map of the spatial distribution of the total integrated 
H I column density in Dw2 is shown in Fig. 4h. The total- 
H I distribution is rather irregular and thus poorly suited 
to estimating the inclination of Dw2. 

The projected velocity field shown in Fig. 4g was deter- 
mined by fitting single Gaussians to each velocity profile 
through the data cube, in a procedure identical to that 
followed for the Dwl data. The observed velocity field 



Table 2. Rotation curves (deprojected for inclination) for 
Dwingeloo 1 and Dwingeloo 2, derived from least-squares mod- 
elling of concentric, rotating, tilted rings 

Dwingeloo 1 Dwingeloo 2 

pX pX 

(") (tolls ') (--) (tons ') C) 

20 20.9±1.6 135.2±5.5 19.5±2.6 296.7±6.0 

40 42.3±1.1 127.3±2.1 28.4±1.2 297.9±2.4 

60 67.8±1.0 115.4±1.1 33.5±1.3 291.4±1.9 

80 83.8±0.7 114.3±0.6 37.9±1.8 281.8±2.7 

100 90.9±0.6 112.0±0.5 43.5±0.9 274.8±1.3 

120 94.1±0.5 111.3±0.4 48.0±0.9 269.5±1.1 

140 97.9±0.5 113.2±0.4 51.6±1.0 271.4±1.1 

160 101.6±0.4 114.2±0.3 

180 104.9±0.5 113.8±0.4 

200 108.3±0.5 113.0±0.4 

220 110.7±0.4 112.4±0.3 

240 112.4±0.4 111.9±0.3 

260 113.5±0.4 111.6±0.3 

280 113.6±0.4 111.8±0.3 

300 112.7±0.5 112.3±0.3 

320 113.9±0.5 113.0±0.3 

340 115.4±0.4 112.4±0.3 

360 114.2±0.5 111.7±0.3 

380 112.3±0.5 112.0±0.3 

400 110.6±0.6 111.2±0.4 

420 110.2±0.6 112.0±0.4 

440 109.0±1.1 114.9±1.0 



seems to be twisted, suggesting a dominant central bar 
or a warp in the outer regions. Modeling the velocity field 
of Dw2 by tilted concentric rings, assuming circular mo- 
tions and thus ignoring a possible bar potential, shows a 
rather constant inclination angle of 69° but a clear change 
in position angle from 294° in the inner regions to 271° in 
the outer ones. The results of the tilted-ring fit are shown 
in the three panels of Fig. 4f. Note that the kinematic po- 
sition angles do not reflect the position angle of the overall 
H I distribution; this situation favors the idea of a bar po- 
tential over that of a warp. 

The left-hand panel of Fig. 4e shows a position- velocity 
diagram at the position angle of the kinematic major ax;is 
in the inner region of Dw2; the right-hand panel shows 
the position-velocity slice along the kinematic minor axis 
in the inner region. Both diagrams are centered on the 
derived dynamical center of Dw2. The dots represent the 
rotation curve projected on to these slices. 

Fig. 4d shows the deprojected Hi surface density, in 
units of M0 pc~^, plotted separately for the approaching 
and receding sides of Dw2. The integrations were done 
in ellipses prescribed by the results from the tilted-ring 
modelling. Although the central H i surface densities in 
Dw 1 and in Dw2 are quite similar, the densities in Dw2 
fall rapidly with distance from its center. 
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Fig. 4. Global kinematic and spatial properties of the Hi distribution in Dwingeloo 2. Panel a) shows, for comparison with 
the global Hi situation, a greyscale representation of the UKIRT A'-band image kindly provided by R. Tilanus and C. Aspin. 
Panel b) shows the CLEANed continuum map, with contours drawn at levels of —1.2, —0.58, 0.58 (the 2-a level), 1.2, 2.3, 4.6, 
9.3, and 19 K. Panel c) shows the global Hi line profile in mjy, corrected for the sensitivity of the primary beam; the vertical 
arrow shows the systemic velocity of Dw2, 94.0 km s~^, as derived from the velocity-field data. Panel d) shows the deprojected 
Hi surface density in units of Mg pc"'^, as determined in tilted rings separately for the approaching and receding sides of Dw2; 
the full-drawn line corresponds to the average of the surface densities measured on the two sides. The panels labelled e) show 
position- velocity maps of the Hi emission observed, at the indicated position angles, along the major axis (left) and along the 
minor axis (right) of Dw2. Contours in the major-axis position-velocity map correspond to —0.21, 0.21 (the 2-a level), 0.42, 
and 0.85 xlO^° Hi atoms cm ^; in the minor-axis map, the contours correspond to —0.20, 0.20 {2a), 0.40, 0.81, and 1.2 xlO^" 
Hi atoms cm~'^. The dots follow the rotation curve as projected onto this position-velocity maps. In both panels, the horizontal 
dashed line gives the systemic velocity while the vertical one corresponds to the location of the center of rotation. The cross 
in the lower-lefthand corner of each of the position- velocity slices gives the spatial and kinematic resolution of the WSRT data 
on Dw2. Panel g) shows the observed velocity field determined for Dw2 by Gaussian fitting to individual Hi spectra. The 
lighter greyscale and the black contours correspond to the approaching side; the darker greyscale and the white contours, to 
the receding. The heavily-drawn black contour follows the systemic velocity. The interval between contour lines is 15 km s~^, 
measured from the systemic velocity of 94.0 km s~^. Panel i) shows the modelled velocity field as a "spider diagram" derived by 
fitting the observed field with tilted rings in which the H I describes circular rotation; the isovelocity contours in this panel are 
drawn as in g). Panel j) shows the difference between the observed and modelled velocity fields; the contours give the velocity 
residuals at levels of —10, —5, 5, and 10 km s^^. Panel h) shows the distribution of the total Hi column density across Dw2. 
As indicated in the text, the noise level varies across this map. Grcyscalcd pixels in the total-H I map all have a signal-to-noise 
ratio > 1.5; the contours correspond to 1.3 (mean S/N = 3), 2.2, 4.2, 6.3, 8.4, and 13 xlO^° Hi atoms cm~^. The three panels 
in /) show the radial variations of, respectively beginning with the upper of the three, the inclination of tilted rings fit to Dw2, 
the position angle of these rings, and the rotation velocity. The shaded circle in the lower left of the two upper panels and in 
the lowermost left-hand panel indicates the FWHM of the synthesized WSRT beam. 



5. Concluding remarks on the intrinsic properties 
of Dwingeloo 1 and Dwingeloo 2 

Figs. 5 and 6 place Dwl and Dw2 together in their relative 
position and orientation on the sky, with Fig. 5 showing 
the total H i maps for the two galaxies and Fig. 6 their 
velocity fields. 

The assignment of Dwl to the class of barred spiral 
galaxies is most clearly required by the appearance of the 
UKIRT K-hand image obtained on August 10, 1994, by 
R. Tilanus, C. Aspin, and T. Geballe. (This image was 
printed in the lAU General Assembly Newsletter "Sidereal 
Times", #3, August 18, 1994.) The K-band data show 
the unmistakable signature of a bar, oriented along the 
major axis of the Hi distribution as determined from the 
WSRT data, and extending to about 1.5 arcminutes from 
the center of the galaxy. The INT and Wise Observatory 
7-band images confirm this assignment (see Loan et al. 
1995). Both the K- and /-band images show spiral arms 
emanating from the ends of the stellar bar, and wind- 
ing over some 180° in azimuth. The superposition shown 
in Fig. 7 of the WSRT total-H i map overlayed on the I- 
band image shows that the optically-traced structure is 
continued to larger radii in the H I distribution. The resid- 
ual velocity-field map shown in Fig. 2j suggests that the 
spiral structure at larger radii perturbs the kinematics of 
Dwl in a manner not uncommon for such galaxies. The 
plot of H I surface density against distance from the center 
shown in Fig. 2d shows a relative deficiency of H i emission 
from the region of Dwl occupied by the stellar bar; such 
morphology is also not uncommon. 



Dwl displays a quite simple morphology. There is no 

indication of a substantial lopsidedncss in the H i surface 
distribution, nor is there any compelling indication of the 
kinematic or spatial distortions which would be expected 
if the galaxy was severely warped, or if it were undergo- 
ing an interaction with a close-by neighboring system of 
comparable or higher mass. 

Dw2 is located approximately along the extension of 
the major axis of Dwl. Unlike the simplicity pertaining 
for the larger galaxy, the kinematic and spatial proper- 
ties of Dw2 do suggest that this smaller system might be 
distorted by an interaction with the more massive Dwl. 
The isovelocity contours shown in Fig. 6 and in Fig. 4g 
show a twist which might be associated with a gravita- 
tional distortion of the smaller system by the larger one. 
The panels in Fig. 4f showing the variation of inclination 
and of position angle of the fitted tilted rings additionally 
suggest a distortion of Dw2. (The panels in Fig. 2f show- 
ing the same information for Dwl indicate no significant 
variation of inclination or of position angle over distance 
from the center of Dwl.) 

Table 2 gives the rotation curves for Dwingeloo 1 and 
Dwingeloo 2, as derived from least-squares modelling of 
the observed velocity fields by concentric, rotating, tilted 
rings; the values tabulated have been corrected for the 
inclinations of the galaxies. We note that the Hi surface 
density of Dw2 drops below dctectability in the outer re- 
gions before reaching a regime of flat rotation. 

The suggested distortion of Dw2 by Dwl remains ten- 
tative until the linear distance between the two systems 
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Fig. 5. Distribution on the sky of Hi column densities measured in Dwingeloo 1 and Dwingeloo 2. The column densities are 
given by the false-color coding as indicated by the color bar on the right-hand side of the figure. The circle in the lower-lefthand 
corner represents the FWHM of the synthesized WSRT beam. Dw2 is located at an angular separation from Dwl of some 21', 
approximately along an extension of the major axis of Dwl. If the two galaxies are located at the same distance from the Sun, 
and if this distance is 3.8 Mpc (see Loan et al. 1995; cf. Kraan-Korteweg et al. 1994), then the angular separation corresponds 
to a linear distance of 24 kpc. 



Fig. 6. Velocity fields revealed by Hi observations of Dwingeloo 1 and Dwingeloo 2. The velocities are given by the false-color 
coding as indicated by the color bar on the right-hand side of the figure. The circle in the lower-lefthand corner represents the 
FWHM of the synthesized WSRT beam. Dw2 is located near the extended major axis of Dwl. There is some indication of a 
kinematic distortion in the smaller system which tentatively suggests that Dw2 is a companion to Dwl. 



Table 3. Intrinsic properties of Dwingeloo 1 and Dwingeloo 2 derived from the WSRT observations 



Dwingeloo 1 Dwingeloo 2 



dynamical center (eq. coord.) 


Q;i960 




02''53'"01".6 


02"50'"21".7 




<5l9BO 




58'*42'43" 


58'*48'07" 


dynamical center (gal. coord.) 


I 




138? 52 


138? 17 




b 




-0?11 


-0?19 


integrated HI flux 




(Jy km s ^) 


205.2 ± 1.2 


30.8 ± 0.5 


HI profile widths, 20% 




(km s-i) 


201.2 ± 0.4 


116.6 ± 1.6 


50% 




(km s-i) 


187.6 ± 0.6 


99.9 ± 1.2 


systemic velocity (LSR) 




(km s-i) 


110.3 ± 0.4 


94.0 ± 1.5 


diameter of HI disk 




(arcmin) 


13.7 


6.4 


continuum flux at A21cm 




(mJy) 


48.2 ± 7.5 


< 16^ 


inclination 




(degrees) 


51 ± 2 


69 ± 3 



3-cr limit for extended continuum emission from Dw2 



is better known. The matter of the distance to the two 

galaxies is, of course, important in several regards. It is 
not reasonable to use the Hubble-flow to gauge the dis- 
tance to such nearby systems. Instead, Kraan-Korteweg 
et al. (1994) estimated the distance to Dwl at about 3 
Mpc on the basis of the blue TuUy-Fisher relationship 
between the extinction-corrected absolute magnitude and 
the inchnation-corrcctcd H i Uncwidth. The galactic ex- 
tinction was determined by Kraan-Korteweg et al. follow- 
ing the precepts of Burstein & Heiles (1982) and using 
the H I column density determined from the data of the 
Leiden/Dwingeloo survey of Hartmann & Burton (1995). 
The principal uncertainties in this estimate stem from the 
large scatter in the blue TuUy-Fisher relationship and in 
the unprecedented use of the Burstein-Heiles relationship 
at the high values of Nm found near to 6 = 0°. It is not at 
all ruled out that the Galactic H i spectra so close to the 
equator are partially saturated (see Burton, 1993), and 
if saturation does cause the extinction derived from the 
integrated H i to be underestimated, then the distance of 
3 Mpc may also be underestimated. Loan et al. (1995) 
derived the distance to Dwl using the /- and i?-band ver- 
sions of the TuUy Fisher relation; this has the advantage 
that corrections for Galactic and internal extinctions arc 
smaller in the near infrared than in the blue. Loan ct al. 
choose as the most likely distance for Dwl the value 300 



km s~^, although they stress that the distance uncertain- 
ties remain substantial. Assuming that the value of the 
Hubble constant is 80 km s— 1 Mpc~^, we are led to adopt 
the value of 3.8 Mpc as the distance to Dwl in order to 
proceed with estimating characteristic physical aspects of 
the galaxies. 

If Dwl is at a distance of 3.8 Mpc then its Hi angular 
diameter corresponds to a linear size of 15 kpc. This es- 
timate of the H I diameter may be compared with the di- 
mension of the ii'-band bar as projected at 3.8 Mpc dis- 
tance, namely about 2.5 kpc, and with the dimension of 
the IRAS 60-/im image (see Loan et al. 1995), namely 
about 4 kpc. The flat regime of the H l rotation curve per- 
tains beyond about 4 kpc. If Dw2 is at the same distance 
as Dwl, its Hi diameter is 7 kpc. The projected angular 
distance between Dwl and Dw2 corresponds to a lower- 
limit on their linear separation of only 24 kpc. 

The total H i mass harbored by the galaxies follows 
from 2.36 x IO^D^/ri, where D is the distance in Mpc 
and Ihi is the integrated Hi flux. For Dwl, Mhi = 7.0 x 
IO^Mq; for Dw2, Mm = 1.1 x IO^Mq. 

The above values for the H i and near- and far-infrared- 
imagc diameters as well as for the H i mass of Dwl are 
consistent with what would characterize a modest barred 
spiral galaxy; the values for Dw2 also seem plausible for 
the small companion. The total mass (i.e. including stel- 
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Fig. 7. Superposition of the WSRT map of the Dwl H I surface densities, shown by contours, on the /-band false-color image 

from the INT (Loan et al. 1995). The Hi surface density contours are drawn at the same levels as those in Fig. 2b; regions of 
highest density are encircled by the red contour. The circle in the lower-lefthand corner represents the FWHM of the synthesized 
WSRT beam. 



lar as well as any dark matter) within a radius r can be 
estimated from Mtot{< f) = {vrotY'i'/G, where Vrot is the 
deprojected rotation speed. For Dwl, Vrot is about 113 km 
s,~^ . Taking r as the distance of the last measured point 
on the rotation curve (see Table 2), namely at about 7f5 
for Dwl or at about 2. '5 for Dw2, then the total mass to 
this distance in Dwl follows at the value 3.1 x 10"'^'' M0 
for Dwl and at the value 2.3 x 10^ M© for Dw2. The ratio 
Mm/Mtot would be about 0.02 for Dwl and about 0.05 
for Dw2; these values seem plausible for a spiral and for a 
gas-rich dwarf, respectively. 

It is interesting to compare the total mass of Dwl with 
that of the Milky Way as well as with that of Dw2; this 
comparison follows simply from the squared ratio of ro- 



tation speeds in the outer flat-rotation regimes, which is 
Mmw/MdwI = (220/113)2 = 3.8 for the ratio of total 
mass of the Milky Way compared to Dwl. For Dw2, where 
H I has not been detected from a flat-rotation regime, we 
take Vrot at the maximum measured value of about 52 km 
s~^. This value then yields a rough estimate of the ratio of 
Mdwi/Mdw2 ^ (113/52)2 = 4.7 for the Dwl/Dw2 com- 
parison of total masses. These ratios give reasonable esti- 
mates of the relative masses within the relevant radii, but 
do not give accurate estimates of the relative total masses, 
because the true total extents are not yet known. At the 
presumed distance and with the total mass estimated as 
indicated, the mass of Dwl would be comparable to that 
of M33. 



W.B. Burton et al.: Neutral hydrogen in the nearby galaxies Dwingeloo 1 and Dwingeloo 2 



17 



Fig. 8. Superposition of the WSRT map of the Dw2 H I surface densities, shown by contours, on the JsT-band false-color image 
from the UKIRT IRCAM3 image kindly provided by R. Tilanus and C. Aspin. The Hi surface density contours are drawn at 
the same levels as those in Fig. 4b; regions of highest density are encircled by the red contour. The circle in the lower-lefthand 
corner represents the FWHM of the synthesized WSRT beam. 



The distance estimate of 3.8 Mpc and the angular 
proximity of Dwl & 2 to the galaxies MafFei 1 & 2 and 
IC342, make the conclusion plausible that these galaxies 
arc all members of the same grouping. Recently McCall & 
Buta (1995) discovered two dwarf companions to Maffei 
1; although these dwarfs would probably contribute little 
to the overall mass budget of the group, their discovery 
additionally motivates searches for galaxies in the vicin- 
ity. We note that before the detection of Dwl there were 
several predictions, for different reasons, concerning the 
likelihood of an undetected galaxy near Maffei 1 and 2. 
Hurt et al. (1993) pointed out a number of morphological 
distortions and other features evident in a deep K-hand 
image of Maffei 2 that were taken as indications of a recent 



interaction with a companion galaxy. Peebles (1990) con- 
cluded that the motions of the Maffei group suggested the 
presence of a nearby galaxy, as then undetected. Whether 
or not the two Dwingeloo galaxies exercise any significant 
dynamical influence on the motions of the Maffei/lC342 
group remains a matter for further consideration. 
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